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Abstract 

The configurational states as introduced by Goldstein represent the system's basins and are characterized by their Iree 
energies ^p{T, V) at a given temperature T and volume V, as we show here. We find that the energies of some of the special 
points (termed basin identifiers here) like the basin minima, maxima, lowest energy barriers, etc., which do not appear in the 
partition function so that the latter is independent of them, cannot be used to characterize the configurational states of the 
system in all cases due to their possible non-monotonic behavior as we explain. The complexity S{^p,T,V), the natural log 
of the number of basins having the same free energy (/^(r, V), represents the configurational state entropy. We prove that the 
equilibrium entropy S{T,V) = S{ipb,T,V)+Sh{T,V), where Sh, and tph are the equilibrium basin entropy and free energy, 
respectively. We further prove that all basins at equilibrium have the same equilibrium basin energy E{T, V) and entropy 
Sb{T, V). Here, ip and E are measured with respect to the zero of the potential energy. The equilibrium complexity is shown to 
be purely an entropic quantity, and not a free energy. The Boltzmann equiprobability principle is shown to apply to the basins 
in that each equilibrium basin has an equal probability V =exp(— 5) to be explored. This principle allow us to interpret the 
basins and their free energy ip as the analogs of the microstates and their energy in the microcanonical ensemble. In addition, 
it allows us to draw some useful conclusions about the time-dependence in the system. We discuss the percolation due to basin 
connectivity and its relevance for the dynamic transition. Our analysis validates modified Goldstein's conjectures that (i) the 
basin free energy measured with respect to the potential energy of the basin minimum is insensitive to the basin minimum 
energy, and (ii) the equilibrium basin free energy iph is insensitive to basins being explored. Thus, our approach demonstarates 
that the configurational state description is not only formally exact, but also analogous to the the conventional approach using 
microstates. All the above results are shown to be valid at all temperatures, and not just low temperatures as originally propsed 
by Goldstein. 



I. INTRODUCTION 

In classical statistical mechanics, the Hamiltonian 
7Y(p,q) = K{p) + E{q) of a system of TV particles in a 
fixed volume y is a sum of the kinetic energy K(p) and 
the potential energy i?(q) of the particles; here p and 
q represent the collective momenta and positions of the 
particles, respectively. The dimensionless total canonical 
partition function (PF) Zn (T, V) of the system can be 
written as a product of two independent integrals 



ZNiT,V) 



1 



(27rfi) 



3N 
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dp 



dq. 



(1) 

Here dp,dq represent integrations with respect to the col- 
lective momenta and positions of the particles, and f3 the 
inverse temperature 1/T in the units of the Boltzmann 
constant /cb- 



A. Kinetic PF 

The first integral Zpf,KE{T) (including the prefactor) in 
the square brackets in (1) contains no interaction energy, 
and is trivial to evaluate because of the quadratic nature 



of K{p). It is clearly independent of the volume V. The 
corresponding free energy Fi\[^ke{T) = —Tin Zn^ke{T) 
is the free energy due to the kinetic energy and is the 
same for all systems, regardless of their potential energy 
of interaction, and volume. The corresponding entropy 
Sn^T) EE ~{dFN,KE{T)/dT) and energy En.MT) = 
Fn,ke{T) + TSn,ke{T) are independent of V, and have 
the same values for all systems. Therefore, we do not have 
to specifically take into account the contribution from the 
kinetic energy part. All we need to do in this case is to 
add the contribution Fn,ke{T) to the free energy from 
the second integral in (1) at the end of the calculation. 
This is precisely what we do here and mostly consider 
the configurational PF discussed below. 



B. Configurational PF 

The dependence on the potential energy (PE) in the 
PF Zn{T,V) is through the second integral Zn{T,V) 
(the prime indicating integration over distinct configura- 
tions of the particles), conventionally called the configu- 
rational PF: 



ZNiT, V) 



dqe 
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(2) 
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The configurations of the system are represented by 
phase points in the phase space TnCV) defined by the 
position coordinates q. Each phase point represents a 



1 



microstate of the system. By appending another axis 
to r7v(^) to represent the potential energy E{q) of 
each point q, we introduce a 3A^-dimensional energy sur- 
face T,pf{V) in a (3A^ + l)-dimensional landscape space 
^n{V). The height of the PE surface represents the po- 
tential energy E of each phase point in the space F n (V) . 
Thus, each configuration or microstate is characterized 
by its potential energy E, along with V, and N; the parti- 
cle number is kept fixed in this work. For the moment, 
let us keep V fixed and consider the slice dSjv(-E', V) of 
the surface with energies between E and E+dE. 

The slice is given by the intersection of 'E,n{V) with 
the two parallel constant energy surfaces T,]\;{E) and 
Ejv(£' + dE) of energy E and E + dE, respectively, in 
A7v(y). This slice defines the microcanonical ensemble 
corresponding to energy E and, according to Boltzmann 
[1], all configurations or microstate are equiprobable to be 
explored by the system. In the canonical, i.e., the NTV- 
ensemble, the entire surface Y^NiV) uniquely determines 
the thermodynamics of the system, and is therefore of 
central importance. This point has been made forcefully 
in the context of glasses in a recent work [2] to which we 
direct the reader. 

In equilibrium, only configurations that have their en- 
ergies E equal to the equilibrium energy E{T, V) deter- 
mine the thermodynamics. These are what we will call 
the equilibrium configurations at T. The equilibration 
process, therefore, is highly selective and picks out only 
the equilibrium configurations that change from one tem- 
perature to another; we are of course neglecting the ther- 
mal fluctuations. It is well known [1] that the energy 
fluctuations AE at constant V,N is given by Ty/Cv,Ni 
where Cv.n is heat capacity at constant V, N. Thus, the 
relative fluctuation AE/E{T,V) behaves as '^/y/N for 
any finite T. This justifies our neglecting the fluctuations 
for a macroscopic system. 

In the following, we will only consider a macroscopic 
system for which N, and V are very large. The ther- 
modynamic limit, which will be implicit in the following, 
requires considering N oo.V —f oo such that n = N /V 
is kept fixed in the canonical, i.e., the A^T^-ensemble. In 
the related A/'TP-ensemble, we only require N ^ oo. We 
will usually suppress the index N (but continue to show 
V) in the following for simplicity. 

C. Goldstein's Ideas and Bcisin Identifiers 

We follow Goldstein and consider the 7VTT/-ensemble. 
The PE surface Tj{V) contains many local minima whose 
importance has been argued by Goldstein [3, 4] for super- 
cooled liquids. He also provided a qualitative description 
of the nature of the resulting landscape created by the 
distribution of minima and drew attention to the idea 
of configurational states, which are associated with PE 
minima in The potential wells containing the min- 

ima are known as basins or valleys. The configurational 
states are collections of allowed configurations within a 



basin. Therefore, the equilibrium configurational states 
must contain only equilibrium configurations of energy 
E = E{T, V). However, because of the importance of the 
basin minima (BM) at low temperatures, Goldstein has 
postulated that the configurational states, instead of be- 
ing associated with the entire basins, are characterized 
by the basin minima (BM), and in particular by their 
energy Ej^. The BM are not the only point of interest 
in basins; the latter have other points of interest such 
as the lowest energy barrier (LB) or the most probable 
energy barriers (MPB), i.e. the barriers with the highest 
population, which are bound to play an important role 
in crossover dynamics or diffusion, as the temperature is 
raised; this point is discussed further in the last section. 
Thus, one can also take the energy E^^ corresponding 
LB or E^^^ corresponding to MPB to characterize the 
basin depending on the temperature range. One can also 
take the highest energy E^ of the basin, the highest bar- 
rier (HB) energy E^^ , or the energy of some other special 
basin point to characterize the basin. In the following, 
we will use the term basin identifier to refer to any of the 
special points of the basins and denote the corresponding 
energy by S. 

It is obvious that each basin identifier has its useful- 
ness limited to a certain temperature range. Thus, as 
the temperature is raised, different identifiers and their 
energies become relevant and to focus on only just one of 
the identifiers like BM will have a limited applicability. 
To provide a description of the system that is useful at 
all temperatures, we need to focus on the basin free en- 
ergy which is determined by the shape and topology of 
the basins. This will require a methodology to express 
the basin free energy in terms of the equilibrium energies 
of an identifier in its temperature range. One of our aims 
here would be to describe such a methodology. 

The shape and topology uniquely determine the free 
energy of the basins at all temperature, which then 
uniquely determine the thermodynamics of the viscous 
fluid. Because of this, we will describe the configura- 
tional states by specifying the basin free energy so that 
the characterization remains valid at all temperatures. 
Later on, we will sec that wc can also use the basin en- 
ergy to represent the configurational states, because the 
equilibrium configurations are identified by the equilib- 
rium energies, as said above. 

The main motivation behind the introduction of con- 
figurational states appears to be the ability to express 
the entropy as a sum of two different contributions at 
low temperatures where the system is confined to one of 
the basins [3]. This partitioning requires the two contri- 
butions to be independent. Could such a partition and 
the independence of its two parts be valid at all tem- 
peratures? In his analysis, Goldstein who only considers 
BM as the basin identifier to characterize the basin has 
listed several conjectures that were common in the field 
[4] at the time and are expected to be valid at low enough 
temperatures. 

(CI) The basin free energy is independent of the possible 



2 



basin being explored. 

(C2) The basin free energy is particularly independent of 
the basin's potential energy minimum Ej^. 

(C3) The partition function is a product of the basin 
and the BM partition functions Zh{T,V), and 
Zbm{T,V), respectively. 

The BM PF Zbm [4] is defined as 

Zbm{T, V)= [ Nbm{El, V)e-('^^dEi^. (3) 



Here, NBM{E-LTV)dEi^ represents the number of poten- 
tial energy minima in the energy range -El and Ej^+dEi^. 
However, the basin PF Zb is left undefined. From the 
form of Zbm in (3), it is most certain that for a given 
basin is defined by considering energies measured with 
respect to the BM Ei^ of that basin; see also Schulz 
[5]. Goldstein emphasizes that basin anharmonicity or 
the curvature at its minimum [6] may be very impor- 
tant and should be included in Zb- According to (CI), 
all equilibrium basins have the same basin free energy 
fb{T,V) = -Tin ^b. According to (C2), fbiT,V) is in- 
dependent of PE minima of the equilibrium basins. This 
conjecture merely reflects the fact that fh{T, V) is mea- 
sured with respect to the basin energy minimum [7] . Ac- 
cording to (C3), the PF [without Zke{T)] is given as a 
product [4, 5] 



Z{T,V) = Zk{T,V)ZBM{T,V). 



(4) 



The form of Zbm in (3) implies the equiprobability hy- 
pothesis that basins having their BM at the same are 
equiprobable. Can these conjectures be justified at all 
temperatures? 

Our aim in this work is to provide a unified formalism 
that remains valid at any temperature. As we will see, 
this will require modifying the above conjectures slightly, 
which can then be shown to be valid at all temperatures. 



1. Condition of Equilibrium and T-Dependence 

Remark 1 For a macroscopic system in equilibrium, the 
PF must be maximized. 

We apply this condition to Zbm{T, V). In equilibrium 
at a given T,V, it should be dominated by basins having 
their minima at the equilibrium minima energy £'l (T, V) 
at which the integrand in (3) is maximum. In terms of 
the BM entropy Sbm{Et^, V) = InNBuiEj^, V), Ej^{T, V) 
is determined by 



{dSBM{Ej^,V)/dEi^) 



V,N 



BL=£L(T,y) 



0, (5) 



which is the condition of equilibrium. The equilibrium 

BM entropy Sbm {T, V) in the canonical ensemble is the 
value of the microcanonical entropy ^bm (-El , V) at -El = 



-El(T, V). Even though S'bm(E'l, V) is independent of T, 
the canonical entropy Sbm{T,V) = Sbm[El{T,V),V], 
obtained after the integrand has been maximized to 
achieve equilibrium and -El replaced by its equilibrium 
value -El (T, V), is a function of T. This is a common 
feature in statistical mechanics. 

Remark 2 A microcanonical quantity is independent of 

T but becomes T-dependent in the canonical ensemble 
through its dependence on some equilibrium quantity such 
as the T-dependence of Sbm{El,V) at El = El{T,V) 
above. 



2. Basins at Equilibrium 

Provided (C3) is valid, the basins that contribute at 
equilibrium have the same equilibrium BM energy -El = 
-El(T, V) and are equally probable. The BM free energy 
Fbm{T, V) = -Tin Zbm(T, V) is given by 

Fbm{T, V) = E^{T, V) - TSbm{T, V), 

and the configurational free energy F{T, V) = 
-Tin Z{T, V) is given by the sum 

F{T, V) = /b(T, V) + Fbm{T, V) = 

[eb(T, V) + E^{T, V)] - T[5b(T, V) + 5bm(T, V)], 

where eb(r, 1/) and Sh{T,V) are the equilibrium basin 
energy (measured with respect to the basin minima) and 
entropy, respectively. The first quantity, therefore, rep- 
resents the equilibrium basin energy measured with re- 
spect to the zero potential. Adding to this Fke(T) = 
Eke{T) - T5ke(T) from the kinetic PF, we obtain the 
free energy 

jr(T,y) = -TlnZ{T,V) 

= E,ib{T, V) - T[5vib(T, V) + Sbm{T, V)], 

where £;vib(T, V) ^ Ekb{T) + eb(T, V) + Ei^{T, V), and 
5'vib(T, V) = 5ke(T) + Sb{T, V) represent the equilib- 
rium basin vibrational energy and entropy, respectively. 

Within the approximation implied in (C3) which is ex- 
pected to be valid at low temperatures, we observe that 
the equilibrium entropy S{T) in the above free energy 
J-"(T, V) or F{T, V) is a sum of two parts, one of which 
is the configurational state entropy Sbm{T,V), and the 
other one is S'vib(T, F) or Sh{T,V), respectively. Thus, 
the approximation in (C3) and the entropy additivity are 
equivalent in the above approach. 



D. A Different Approach 

According to Goldstein, we can interpret the thermo- 
dynamics of the system as being confined to any of the 

many basins; this is certainly expected to be valid at 
least at low temperatures. His seminal ideas based on 
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the above conjectures (C1-C3) have been central in un- 
derstanding glasses or amorphous materials, as they pro- 
vides an intuitive rational for the observed relaxation 
of two different time-scales [3, 8]. Using these ideas, 
Goldstein was able to provide an interesting and suffi- 
ciently tractable scheme for calculation discussed above, 
and made some quantitative predictions [3, 4]. Our in- 
terest here is to check the validity of the above conjec- 
tures, to see if they are valid at all temperatures, and 
to shed new light on some of the core concepts, like the 
basin free energy, the nature and significance of configu- 
rational states and their relationship with the basin en- 
ergy or free energy, the d priory equiprobability of the 
basins in equilibrium and their energies and entropies, 
etc. We also discuss the phenomenon of percolation of 
basin connectivity and its possible connection with the 
dynamics in the system. 

We carry out an analysis by taking an approach which 
has appeared in the literature [9-12], but has not been 
studied extensively. In particular, we measure the basin 
free energy with respect to a common reference energy, 
the zero of the potential energy (see below) and denote 
it by Lp\, (T, V) to distinguish it from the previous free 
energy f\,{T, V). The approach leads to the equiprobabil- 
ity theorem in a straightforward manner. This is highly 
pleasing as it brings forth the similarity with the conven- 
tional equiprobability hypothesis due to Boltzmann [1]. 
We substantiate the first two conjectures of Goldstein in 
slightly different forms; the third conjecture also needs 
to be modified. All these are shown to be valid at all 
temperatures. 



1. Complexity 

As said above, Goldstein's ideas seem to justify a well- 
accepted belief in the glass community that the entropy 
of the supercooled liquid is a sum of two contributions: 
(i) the entropy ^bm (T, V) due to different equilibrium 
configurational states or basins, and (ii) the entropy 
Sh{T,V) associated with an equilibrium basin. (The 
supercooled liquid is a metastable liquid, which is de- 
scribed by the landscape in which some of the basins 
that are relevant for the crystal phase are deleted from 
consideration.) However, whether such a partition actu- 
ally occurs in viscous liquids has never been rigorously 
demonstrated. (As we will see here, the actual rigor- 
ous demonstration is not trivial.) Borrowing ideas from 
the spin-glass field, Mezard and Parisi [9] introduced the 
idea of complexity S{T, V) as a candidate for the en- 
tropy of equilibrium configurational states, which has 
been adopted by other workers in the field; see, for ex- 
ample, [11, 12]. The complexity is determined not by the 
EM energies E'l (or the energies of other basin identifiers 
such as E^^,E^^^, etc.) but by the basin free energies; 
see below. Coluzzi ct al [10] gave a plausibility argument 
that the liquid's equilibrium entropy S(T, V) is a sum of 



the complexity S{T, V) and Sh{T, V) : 

S{T, V) = S{T, V) + 5b(T, V). (7) 

Coluzzi and coworkers [11] have given a numerical scheme 
to evaluate complexity directly by making several hy- 
potheses, one of which is the equiprobability of the basins 
of the same free energy, and another one is the above en- 
tropy partition (7). 

2. Present Goal 

Since it appears, from all accounts, that the central 
idea of configurational states or basins and the complex- 
ity are important concepts that play a pivotal role in the 
thermodynamics at low temperatures and hopefully help 
us unravel the mystery of glassy states, it is also impor- 
tant to demonstrate that the above partition in (7) is 
valid rigorously and that the equiprobability hypothesis 
can be justified on general grounds not only at low tem- 
peratures but at all temperatures. The latter allows us to 
treat the basin free energy ifi{T, V) in the same manner 
as the energy in the micro canonical ensemble because of 
the equiprobability of configurations of the same energy 
in that ensemble. The partitioning of S{T, V) replaces 
the third conjecture of Goldstein. 

We will prove the following statements. 

(51) The basin free energy fj{T,V) of the j-th basin 
is independent of the energy of a particular basin 
identifier, like the BM energy. 

(52) The equilibrium basin free energy ipi,(T,V), energy 
Eiy{T.V), and entropy Sh{T,V) are the sam,e in 
all equilibrium basins that are explored at a given 
T,V. 

(53) All equilibrium basins have equal probability of being 
sampled. 

(54) The entropy partition {7} is a general result. 

We discuss the general topological features of the land- 
scape in the next section, and write down the fundamen- 
tal thermodynamic differential relations. We then com- 
pare this analysis with the standard analysis, which is 
carried out in Sect. III. The results are discussed in the 
last section, which also contains some speculative com- 
ments and a short summary of our results. 

II. GOLDSTEIN'S LANDSCAPE 
A. Analysis 

The surface S(F) can be formally decomposed into 
nonoverlapping basins Bj{V). There would be an en- 
ergy barrier between different neighboring minima. In 
each basin Bj{V), the allowed energy range is AjE = 
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EY — Ei,j, where and Ei,j{V) are the highest and 
lowest ahowed energies in Bj{V) [13]. The highest energy 
also represents the highest barrier energy in the basin. 
The exceptions are the basins that contain particles near 
the walls of the container (of volume V) in which the 
interactions of the particles with the walls will make the 
potential very large oo) in order to confine the parti- 
cles. For such basins, Ej oo. In this case, the highest 
energy does not represent the highest barrier to some 
other basin. The lowest energy in a basin does not nec- 
essarily have to be point-like in that it is associated with 
only one phase space point q. It is possible that there are 
many contiguous phase space points, all having the same 
lowest energy, with no energy barrier between them such 
as the minima in an inverted Mexican hat. In this case, 
all these points would be classified as representing a sin- 
gle minimum and belonging to a single basin. The same 
situation may also occur for other points like the points 
of highest energy, the cnca-gy barriers, etc. The number 
of basins in r(V") will be denoted by Nb{V). The canoni- 
cal partition function Z{T, V) can be formally expressed 
in the form of an identity as a sum over the basins: 



Z{T,V)^ J2 Z,iT,V), 

j=l,NB{V) 



where 



Zj{T,V) 



dqe 



-/3E(q) 



(8) 



(9) 



is the PF of the j-th basin Bj{V), and the integration 
is over the coordinates confined within the basin. The 
energy is measured with respect to a common zero of po- 
tential for all basins, and not with respect to some prede- 
termined £j (V) associated with a basin identifier. We de- 
note the basin free energy by fjiT, V) = —TlnZj{T, V) 
and introduce fj{T,V) = ifij{T,V) — 5^(1^), which rep- 
resents the basin free energy measured from £j{V). The 
probability Vj (T, V) that the system under equilibrium 
explores the basin Bj{V) is given by 

V^iT,V)^Zj{T,V)/Z{T,V). (10) 



B. Fundamental Thermodynamic Relations 



1. Basins 

Let us rewrite (9) in a different way as follows. Con- 
sider the surface T,j{V), which is the part of the sur- 
face Ti(V) that belongs to the basin Bj{V). The surface 
I]j (F) exists for energies in the range ^jE. We project 
it on the 37V-dimensional phase space r(V^), which is 
orthogonal to the E-axis. The projection Ilj{V) deter- 
mines the number of configurations (not to be confused 
with the configurational states of Goldstein) of the sys- 
tem that belong to the basin as follows. Consider the 



slice d'Sj{E,V) of the surface 'Sj{V) between the ener- 
gies E and E + dE, E e ^jE. The slice contains all the 
points in Tij{V) that lie between the intersections (known 
as the turning points in classical mechanics) of Tij{V) 
with the two parallel constant energy surfaces S(£^) and 
Y,{E+dE) that were introduced above, and will be closed 
(in the sense that a ring is) if the slice cuts the basin in 
all directions; otherwise it is open. For energies close to 
E-Lj{V), the slice is closed, but becomes open as we get 
to higher energies and eventually disappears at E^IV). 
The projection dIij{E, V) of this slice on T{V) is closed 
(open) if dT,j{E, V) is closed(open). The area \dIlj{E, V)\ 
of the projection determines the number of configura- 
tions Wj{E, V)dE = \dnj{E, V) \ between the energies E 
and E + dE. Let us introduce the microcanonical hasin 
entropy function Sj{E, V) = \nWj{E, V), which is inde- 
pendent of T. The entropy Sj {E, V) at fixed V, N satisfies 
the differential relation 

dSj{E, V) = {dSj{E, V)ldE)y,, dE. (11) 
Now, we can rewrite (9) as follows: 

Zj (T, V)= j dE exp[5j {E, V) - 13E{V)] . 

For a macroscopic system at fixed T, the integral is dom- 
inated by the maximum integrand at £■ = Ej (T, V) (see 
Remark 1), given by 



{dS,iE,V)/dE) 



V,N 



E=Ej(T,V) 



= /?, 



(12) 



which is the standard thermodynamic relation between 
the entropy and temperature; compare with (5). The 
energy Ej (T, V) represents the average basin energy at 
a given temperature T,V, (and N). (We will reserve the 
usage of the term "equilibrium" for a thermodynamic 
quantity only after the PF Z{T, V) has been evaluated. 
To make this distinction, we use the term "average" at 
this stage.) Let us introduce the energy difference AEj = 
E-£j, and A£;j (T, V) = Ej (T, V) - £j for the j-th basin, 
where £j is the energy of some particular basin identifier 
like BM, LB, etc. We can rewrite (12) as 



idSjiE,V)/dAE. 



3>V,N 



AEj=AEj{Ty) 



p. (13) 



The canonical basin entropy Sj (T, V) is given by the 
microcanonical basin entropy evaluated at the average 
E = Ej{T,V): 

Sj{T,V) ^ Sj[Ej{T,V),V]. 

Using (11) &t E = Ej{T, V), and (12), we find that 

dEj (T, V) = TdSj (T, V) (14) 

at fixed V, N. The average basin free energies ipj (T, V) 
and fj (T, V) = (pj (T, V) — £j are given by 



^,{T,V) 
fi{T,V) 



E,{T,V)-TS,{T,V), (15a) 
AE,iT,V)-TSj{T,V), (15b) 
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A given T corresponds to a particular value of AEj (T, V) 
[see (13)] measured from £j. Thus. AEj{T, V) will remain 
the same even if the value of £j is changed by shifting 
the origin of the energy. In other words, AEj {T, V) and 
the average basin free energy fj (T, V) do not depend 
on £j. This proves our statement (SI). Consequently, 
i^j{T,V) = fj{T,V) + £j depends linearly on £j. Using 
(11), (12), and (14), it is easy to obtain the following 
fundamental differential relation at constant V,N : 

d^j (T, V) = dfj (T, V) = -Sj (T, V)dT. (16) 

This gives the change in the basin free energy due to the 
temperature change within the basin Bj{V). 



not necessary in (12) since Sj{E, V) did not have an T- 
dcpcndcncc, even though the variation there was also car- 
ried out at fixed T. The equilibrium complexity S{T, V) 
is given by S{(p, T, V) evaluated at (p = (pb{T, V) : 

S{T,V) = S[MT,V),T,V]. (20) 

(Here, wc have used an overbar to indicate that S{T, V) 
is a function of two variables, while S[ipb(T, V), T, V] is 
a fiuiction of three variables.) Each equilibrium basin 
has the same free energy (^b- This proves a part of our 
statement (S2). 

3. Basin Entropy and Complexity 



2. Complexity Analysis 

At a given T, V (and N), the basin free energy Lpj (T, V) 
varies from basin to basin; hence, (pj (T, V) represents a 
family of functions, one for each j. Let = <p>{T, V) 
represent one of the various functions ipjfs in this family, 
and Af{tp, T, V) the number of basins having the same 
free energy ipi for a given T, V, regardless of whether £j is 
the same or not for these basins. In other words, these 
basins need not all have the same basin identifier energy 



£j = £. We now rewrite (8) as 

Z{T, V) = Y,N{^,T,V)e 



■I3<p 



(17) 



For a macroscopic system, we expect the family of func- 
tions represented by (p{T, V) to be almost continuous 
"over" the basin index j. This is an assumption that is ex- 
pected to remain valid even at low enough temperatures. 
With this assumption, the sum in (17) can be replaced 
by an integration over this family. 

Let us introduce the complexity S{ip,T,V) = 
lnJV{ip, T, V), which satisfies the following differential re- 
lation at constant V (and N) : 

dS{ip, T, V) = {dS/d^)^^y # + {dS/dT)^ y dT. (18) 

For a macroscopic system at fixed T,V, see Remark 1, 
the PF in (17) is dominated by that particular value of 
the basin free energy, which we denote by iph = ^h{T, V), 
for which the summand (or integrand) is maximum, i.e. 
<&(</?, T,y) = — TS is minimum. The basins having 
this free energy </?b {T, V) represent the equilibrium basins 
determining the thermodynamics at fixed T, V. As said 
above, these equilibrium basins are not required to have 
the same basin identifier energy £j = £. The location of 
the maximum at if = (pb, for continuous (p, is determined 

by 



{dS{ip,T,V)/d^).r,v 



(19) 



where the derivative is obtained by a variation over the 

family of function (p. The temperature T is kept fixed un- 
der this variation and is shown explicitly in (19); this was 



We can invert the function S{(f, T, V) over the entire 
positive temperature range where it is monotonically in- 
creasing with (p, sec (19), and express </? as a function 
of S : ip = f{S, T, V) over this range. Accordingly, the 
change in (p due to changes in the complexity <S and T is 
given by 

dip{S, T, V) = {dip/dS)j. y dS+ {dip/dT)^^ dT. (21) 

In equilibrium, ip and S should be replaced by p^ and 5, 
respectively, and the coefficient of dS, according to (19), 
should be T. Prom (18), we note that at equilibrium 

{dMT, V)/dT)sy = -T {dS{T, V)/dT)^^ y . 

In the identity 



di£b 
dT 



s,v 



+ 1^ 
dS 



T,V 



dT 



obtained from (21), the first term on the right represents 

the change in the equilibrium basin free energy, see (16), 
solely due to the temperature change dT within a basin; 
the contribution from the complexity change due to dT 
is given by the second term. Thus, the first term is the 
negative of the basin entropy S'b(T, V), and we have 



dTjy 



T 



dT 



Sb{T, V). 



(22a) 
(22b) 



4- Free Energy and Entropy 



The equihbrium free energy F{T, V) = -TlnZ{T, V) 
is given by the value of ^{<p>, T,V) = (p — TS at equilib- 



rmm, i.e.. 



F{T, V) = $((^b, T, V) = pbiT, V) - TS{T, V). (23) 

Taking its differential at constant V, (and N), we obtain 

dF{T, V) = diphiT, V) - TdS{T, V) - S{T, V)dT 

= -[S{T, V) + 5b(T, V)]dT, (24) 
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where we have used (21) appHed to (^b- The entropy can 
be calculated directly from S{T, V) = -{dF/dT)v From 
(24), we have 



S{T,V)=S{T, y) + 5b(T, V) 



fdTS\ 



where we have used (22a). Thus, the entropy of the sys- 
tem is a sum of two contributions, which proves our state- 
ment (S4). The derivative on the right-hand side in the 
above equation can also be written as — {dF/dT)^^ y , 
which should be identical to the derivative —{dF/dT)v 
from the definition of S. Thus, 



{dF/dT)v = (dF/dT) 



which shows that the equilibrium free energy F{T, V) 
is independent of the basin free energy (/?b (at fixed T). 
This should come as no surprise since the equilibrium ip^, 
is determined by the condition 



{dF/dip^)T,v = (5$(^,r, v)/d^)T,v = 0. 



(25) 



We can now include the kinetic energy part of the free 
energy to obtain the complete free energy 

^(T, V) = -TlnZ(r, V) - {Sb(r, V) + Eke{T) 
- T[Sb{T, V) + Ske{T)]} - T5(T, V). 

The quantity in the curly brackets represents what is con- 
ventionally called the equilibrium vibrational free energy, 
even though one can be sure of its vibrational nature 
only at low temperatures where we expect the projec- 
tions dHj {E, V) of each contributing basin to be closed 
(for vibrations, the basin must confine the motion in all 
directions, which may not be possible at high temper- 
atures). The last term is the contribution due to the 
complexity. 



5. Equilibrium Identifier Energy £(T,V) 

The equilibrium basin identifier energy E{T,V) is, by 
definition, given by 

£{T,V)^ J2 £jiV)VjiT,V). (26) 

j = l,NBiV) 

It is clear that this equilibrium value cannot be expressed 
as a derivative of the partition function. It is easily reme- 
died. We add a dimensionless parameter A that couples 
to the basin identifier energy £j so that there is an ad- 
ditional term X£j in the energy, and define a modified 
PF 

Z{X,T,V)= e-^^^i^^^Zj{T,V). (27) 

j=l,Ns(V) 

The equilibrium £{T, V) is given by 

£{T,V) = -T[{d\nZ{X,T, V)/dX)T,v]x=o- 



In the limit A 0, this modified PF reduces to the earlier 
PF in (8). For A > 0, basins with lower identifier energies 
are weighted more than the higher ones, with the reverse 
true for A < 0. For A = 0, there is no preference for 
any particular value of £j ; recall that Zj {T, V) does not 
depend on £j. It is obvious that the sum in (26) can be 
limited to only equilibrium basins for a macroscopic sys- 
tem. For equilibrium basins, all allowed values of £j are 
equally probable and will result in the maximum spread 
A£ in the possible values of £j under the summation in 
(26). Thus, it appears that the spread in the allowed val- 
ues of £j may very well be of the order of N, and not y/N. 
In other words, the summand in (26) is most probably 
not peaked sharply at some particular value of £; see also 
Sect. IV(D) below. We hope to answer this question in 
future. 

It is easy to show from (26) that 



dT 



Nb{V) 

J2 A£j{V)AEj{T,V)Vj{T,V), (28) 

i=i 



where A£j{V) = £j{V) - £{T,V), AEj{T,V) = 
Ej{T,V) - E{T,V), and E{T,V) is the equilibrium 
energy of the system. Thus, the derivative is a 
cross-correlation between two fluctuations, A£j{V) and 
AEj (T, V) . Since cross-correlations do not usually have 
a fixed sign, there is no theoretical reason for £{T, V) to 
be a monotonic increasing function of T. This does not 
mean that the above derivative cannot be positive for 
many physical systems. 

We now show the expected result that the equilibrium 
free energy F{T, V) cannot depend on the identifier en- 
ergy explicitly by demonstrating that 



{dF/d£)T,v = 0; 



(29) 



compare with (25). For this, we use ^{(fi,T, V) = ip — TS 
and obtain 

{d<^/d£)T,v = {d(p/d£)T,v - T{dS/d£)T,v, 
which, from (21), reduces to 

{d^/d£)T,v = {dip/dS)T,v{dS/d£)T,v - T{dS/d£)T,v. 

In view of (19) at equilibrium, the right hand side is 
identically zero, thereby proving (29). 

6. Monotonic £{T,V) 

Because of the importance of different identifiers in dif- 
ferent temperature ranges, it is desirable to express ther- 
modynamic quantities in terms of £{T, V) in these tem- 
perature ranges. Therefore, the question arises whether 
it is possible to relate thermodynamic quantities with 
£{T,V). However, as discussed above, these identifiers 
do not determine the thermodynamics developed here; 
see (29), because of which the relationship with £{T, V) 
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will not be unique, unless some particular scheme is used. 
We present a way to accomplish this under the condition 
that £{T, V) is a monotonic function of T in its temper- 
ature range. We will, therefore, consider this case below. 
The monotonic function £{T, V) of T can be inverted to 
yield 

T = T[£{T, V),V]. 

This will allow us to express any function of T as a func- 
tion of £{T^V). For example, we can express fj{T,V) 
as a function fj[£{T, V), V] of £{T, V). We can similarly 
express ipj (T, V) as a function (pj [£{T, V), V] of £{T, V). 
We can also express (p]^{T,V), and S{T,V) as a func- 
tion <fi,[£{T,V),V], and S[£{t,V),V], respectively, of 
£{T,V). 

Let N-Bi{£ ,V)d£ represents the number of basins in 
which the basin identifier has its energy in the range £ 
and £+d£, and introduce the corresponding microcanon- 
ical entropy Sbi{£, V) = lniVBi(f , V). The T-dependent 
basin identifier entropy evaluated at £ = £{T, V) 

Sbi{T,V) = Sbi[£{T,V),V] (30) 

does not necessarily represent the entropy due to equi- 
librium basins for which £ = £{T,V). The equilibrium 
basins satisfy the conditions in (S2), but the basins enu- 
merated by Sbi{T,V) have no such conditions imposed 
on them. Moreover, the equilibrium basins have a wide 
range of allowed basin identifier energies £j because of 
A = 0, as discussed above. Thus, the entropy due to equi- 
librium basins with £ = £{T, V) is most probably very 
different from the complexity S[£{T, V), V] obtained by 
the above prescription from S{T, V). It is clear that while 
the basins counted by ArBi(T, V) = Nbi[£{T, V),V] con- 
tain all the equilibrium basins with £ = £{T,V), they 
also contain many other basins not satisfying (S2). Let 
^Bi,NE(r, V) represent the latter non-equilibrium basins. 
Since not all equilibrium basins in 77{T, V) have £ = 
£{T,V), the number of equilibrium basins Nbi,e{T,V) 
in Nbi{T,V) is such that Nbi,e{T,V) < JI{T,V). By 
definition, we have 

Nbi{T, V) = 7Vbi,e(T, V) + NbiMT, V). (31) 

The current analysis cannot answer whether the sec- 
ond contribution is thermodynamically insignificant com- 
pared to the first contribution in (31); see the discussion 
of (40) below. To have Nbi{T, V) = J7{T, V) at all tem- 
peratures would be a precarious coincidence and hard to 
justify. Therefore, the sum 

5Bl(r, V)+Sy-,{T, V) 

will most probably not represent the equilibrium entropy 
S{T,V) at all temperatures. For £ =Ei^, Sbi{T,V) rep- 
resents the entropy Sbm{T,V) introduced above in (3). 
Thus, most probably 

S{T, V) ^ Sbm{T, V) + 5b(T, V) (32) 



at all temperatures. In other words, Sbi{T,V) is most 
probably not equal to S{T, V). 

As said above, the representation a thermodynamic 
quantity in terms of £{T,V) is not necessarily unique. 
Consider, for example, S{T,V). There is an alterna- 
tive way to express it in terms of £(T, V) by recalling 
its fundamental definition 'S{T,V) = S[(pb{T,V),T,V] 
given in (20). We express <fh{T, V) in the form (ph{£, V) 
to obtain an alternative expression S^^{£,T,V) = 
S[ipi,{£,V),T,V] for the complexity family in which 
£{T, V) has been represented simply by £. Let us now 
consider the following derivative 

{dS/d£)T,v = {dS^^{£,T,V)/d£)T,v = (3{d^b/d£)T.v, 

(33) 

which can be obtained from (18), and where we have also 
used (19). We should recall that £ = £{T,V) above. 

7. Equiprobability Principle 

Prom (10), and (23), we note that 

Vj{T,V) = exp[-S{T,V)], (34) 

for equilibrium basins, so that the system has equal prob- 
ability exp[— »S(T, V)] to belong to any of the equilib- 
rium basins in Af{(p, T, V). This proves (S3). The basins 
counted in N{(ph , T, V) may have different fj (T, V) = 
(fib{T,V)— £j, as there is no guarantee that all these 
basins will have their basin identifier at the same en- 
ergy. The probability Vj (T, V) is determined only by 
fbiT, V), the free energy of each of the equilibrium basins 
in (23), which in turn determines F{T,V) and S{T,V). 
Comparing with the standard Boltzmann equiprobability 
principle, we can understand (34) by the following anal- 
ogy (see also [11]): the cquilibrimn basins and their free 
energies play the role of microstates and their energies 
in the microcanonical ensemble; the basins and their free 
energies play the role of microstates and their energies in 
the canonical ensemble. 



III. STANDARD PICTURE 

We can rewrite the configurational PF in (2) slightly 
differently in a form, which is more standard and does 
not involve individual basins: 

Z{T,V) = J dEW{E,V)e->^^^^\ (35) 

where W{E, V)dE represents the number of microstates 
or configurations represented by the the slice dT,{E, V) 
with energy between E and E + dE. We obtain the num- 
ber of microstates in the slice dT,{V) by considering the 
area of its projection dn{E, V) on T{V): W{E, V)dE = 
\dn.{E, V)\.lt is evident that 

dU{E, V) = UjdUj{E, V), (36) 
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where the union is over only those basins that do not 
yield null projections; these are the basins for which E ^ 
AjE. The tilde over the union denotes this restriction. 
Consequently, we have 

W{E,V)=^^^^^^^^^^W,iE,V), (37) 

where the tilde has the same meaning as in (36). In- 
troducing the microcanonical entropy S{E, V) = In 
W{E,V), we conclude (see Remark 1) that the PF in 
(35) is dominated by the equilibrium energy E{T, V) for 
which the integrand is maximum, and gives the equilib- 
rium state. The equilibrium energy is the solution of 

[dS{E,V)/dE]^^^^^y^=P, (38) 

which is obtained by the variation of the integrand at 
constant T, V. At a given temperature T, only those 
microstates with energy in the range E{T,V), and 
E{T, V)+dE determine the thermodynamics through the 
entropy S{T, V) = S[E{T, V),V]. The canonical entropy 
S{T, V) along with E{T, V) determine the equilibrium 
free energy F{T, V) = E{T, V) - TS{T, V). 

As the kinetic energy is not a part of the energy E, 
there is an upper bound to the energy [13]. Therefore, the 
microcanonical entropy S{E, V) near this upper bound 
will be a decreasing function of E. This region corre- 
sponds negative temperatures according to (38). Since 
S{E, V) is an increasing function near the lower end of 
the energy E, there is expected to be a maximum in 
S{E, V) as a function of E. The rising branch of S{E, V) 
corresponds to positive temperatures and it is this branch 
that is of interest in this work. 



A. Connection with Landscape: Selectivity 

Using (15a) for the basin free energy (ph{T, V) in (23), 
and comparing with the above definition of the free en- 
ergy in the standard picture, we find that 

E{T,V) = E^{T,V). (39) 

The significance of this result is the following. All config- 
urations contributing to the canonical entropy S(T, V) in 
the thermodynamic limit have the same equilibrium en- 
ergy equal to E{T,V). The above equation (39) shows 
that these configurations, which can also be classified as 
belonging to the different basins, as shown in (37), not 
only have the same equilibrium basin energy Eh{T,V), 
but also have the same basin free energy (p\y{T,V); see 
(20). From (15a), this also implies that they all have the 
same basin entropy Sh(T, V). Thus, the landscape ther- 
modynamics is highly selective. At a given T, only those 
basins contribute to the thermodynamics that have the 
same Eh{T, V), Sb{T, V), and ipb{T, V). This proves our 
statement (S2) in totality. 



B. Equilibrium Basins 

Let A''b (T, V) denote the number of basins for which 
E — i?b(r, F) e Aj i?. We say that theses basins exist 
at T. Not all these basins will satisfy all the conditions 
in (S2). The equilibrium basins are the only basins that 
determine the thermodynamics and must satisfy (S2). 
The number of these basins is A/°(r, V) = N'{ipb,T, V). 
Thus, we expect Jj'{T,V) < Nb{T,V). The munber of 
non-equilibrium basins Ni^e{T, V) that do not contribute 
to the thermodynamics of the system but exist otherwise 
in the above sense is given by their difference: 

Nb{T, V) = W{T, V) + N^e{T, V). (40) 

Based on general arguments, there is no way to answer 
whether the second contribution is tliermodynamically 
insignificant compared to the first contribution in (40). 
All we can be sure of is the following. Because of (7), we 
know for sure that the equilibrium number of configu- 
rations W{T,V) = W[E{T,V),V] = exp[S{T,V)] is ex- 
haustively coming from configurations in all the A7'(T, V) 
equilibrium basins: 

W{T, V) = Jf{T, V) exp[5b(T, V)]. 

The number of configurations coming from the 
A^NE(r, V) non-equilibrium basins is thcrmodynamically 
insignificant. This result has no bearing on the relative 
magnitude of the two terms in (40). It is possible to 
have iVNE(r, V) significantly larger than Jf{T, V) in (40) 
and still have the contribution from the non-equilibrium 
basins to W{T, V) insignificant compared to W{T, V). 

IV. DISCUSSION AND CONCLUSIONS 
A. Basin Connectivity and Percolation 

We have already noted, see (36), that the projection 
dIi{E, V) is the union of different non-null components 
dI\j{E, V). A component may have more than one dis- 
joint pieces, such as the one from a basin in the form of 
an inverted Mexican hat discussed earlier. Correspond- 
ing to some energy E is the equilibrium temperature T 
such that 

E{T, V) = E. (41) 

Because the heat capacity Cy is non-negative, E{T, V) 
is a monotonic increasing function of T over the posi- 
tive temperature range, which is what we are interested 
in here. Thus, increasing (decreasing) E is the same as 
increasing (decreasing) T. 

The number of components at E in the union in (36) 
is Nb{E,V) introduced earlier. The shape and size of 
the components and their number change with E. The 
components dIlj(E,V) cannot be open; they arc either 
connected together to form a cluster, or are isolated, and 
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must be closed (in the form of a ring). Two components 
are connected, when there is no energy barrier separat- 
ing them [15]. As we reduce E or T, these components 
become disconnected due to the presence of an energy 
barrier separating them. As Goldstein has pointed out 
[3], these energy barriers become very important at lower 
temperatures as we will discuss further in the following. 
If a component is isolated, this will mean that there is an 
energy barrier separating this basin from all other basins. 
At very high energies {E > Eqp; see [13]) that are higher 
than all energy barriers, all projections are open, but 
connected together. It contains only a small number of 
basins containing particles near the walls. These energies 
are not of any interest here. As the energy is reduced, 
some energy barriers are cut by the surfaces E(i?) and 
T,{E + dE), which produce isolated closed components 
and will increase the number of components. In the same 
process, some of the components disappear as their basin 
minima energy Ei^ become higher than E. This will re- 
duce the number of components. Thns, we expect the 
number of components to change with E. 

For a macroscopic system, the above idea of basin con- 
nectivity in the 3A/'-dimensional phase space must be 
phrased in a way that makes statistical sense. A basin i 
is connected to a basin j, if the measure of their common 
opening (along which there is absence of any barrier) is a 
non-zero fraction of the projection measure \dni{E, V) \ . 
We will call such openings relevant. The non-zero frac- 
tion of the projection measure will add a non-zero en- 
tropic contribution to the barrier height. (It is almost 
impossible for the system to find a barrier along one of 
the directions out of macroscopically large 37V directions 
over any finite, no matter how long, time; the entropic 
contribution for such a barrier is zero.) Thus, it is the 
free energy , and not the energy that becomes relevant 
for barrier hopping. A relevant opening will have a non- 
zero probability Pij of a jump Bi Bj from the i-th 
basin to the j-th basin. Openings that are zero fraction 
of the projection will be called irrelevant as Vij =0. It 
is in the sense of a relevant opening that the two basins 
will be considered connected below. Connected basins al- 
low the system to probe them without encountering any 
energy barriers. In the following, we will make distinc- 
tions between two kinds of barriers [15]. Barriers whose 
heights above E{T) are of the order of T can be over- 
come due to thermal fluctuations (AE = T^jCy) and 
the time required for the jump Bi Bj is less than or 
equal to the microscopic relaxation time. If the interest 
is to study dynamics at a much longer range of time, we 
can treat these barriers as not different from an opening. 
The barriers whose heights are much higher than T will 
give rise to activation processes requiring much longer 
times to jump over the barriers. These high barriers are 
the ones wc consider in the following. The most probable 
barriers are those barriers that have the maximum prob- 
ability to be overcome at a given T. The barriers whose 
heights above E{T) are of the order of T will be treated 
as openings; thus, they are suppressed [15]. 



It is obvious that the physics of a continuum model 
requires the system in the diffusive regime to probe vari- 
ous configurations that continuously transform into each 
other, as long as we are far away from any phase transi- 
tion or we are in a metastable state obtained by avoiding 
a transition like the crystallization. Therefore, we ex- 
pect various basins to be continuously connected to each 
other in the sense that there is at least one continuous 
passage from one basin to another through intermedi- 
ate basins, without any energy barriers to overcome, so 
that the phase point can move continuously from one 
basin to the other. This will certainly be the case at 
high enough energies (that correspond to high enough 
temperatures). If it were not the case, the system will 
have to overcome energy barriers using activation pro- 
cesses because of which, it will "jump" from the equi- 
librium configuration in one basin into another equilib- 
rium configuration of the other basin (with energy bar- 
riers between them); the configuration of the new basin 
is different from that in the old basin in that the new 
configuration is not in the immediate vicinity of the old 
one. The high-temperature situation will correspond to 
all the components in the projection dIl{E, V) being con- 
nected together as one connected cluster. We can visual- 
ize this situation as the percolation of the basin slice pro- 
jections dIlj{E, V) in one macroscopic cluster, called the 
percolating cluster V{E, V) that spans the entire space 
r(F). Here, the percolation occurs between components 
dIij{E^ F), all of different shapes and sizes. 

The entropy of the system S{E, V) is determined by 
the percolating cluster ViE, V) at high enough energies. 
There are other finite clusters or isolated components in 
addition to I?(i?, V). But not all components and clusters 
are thermodynamically significant and contribute to the 
entropy. Only those components contribute that have 
the same E{T,V),S\,{T,V), and ipb{T,V). Because of 
the continuity in the diffusive regime, we expect these 
components to be part of V^E, V). The number of such 
components is given by N'{iph,T,V). At high energies, 
'D{E, V) contains disordered configurations pertinent for 
the liquid state. The percolating cluster allows the sys- 
tem to probe all of the connected equilibrium basins that 
occur in the diffusive regime. 

As we reduce the energy, some of the basin compo- 
nents (i) disappear since they no longer exist {E < Ei), 
(ii) appear since their energy range contains the energy 
E {Ei^ < E < E^), or (ui) begin to get disconnected 
from others due to the emergence of high energy barriers, 
but the percolating cluster 'D{E, V) persists in dn{E, V) 
although it is of a different size. It may happen that 
some basin components get completely isolated from each 
other and the percolating cluster. It is highly plausible 
that in the interesting range of energies, 'D{E, V) will de- 
crease in size as E is reduced. However, the entropy of the 
system is still determined by the percolating cluster, at 
least at higher energies that correspond to temperatures 
much higher than the melting temperature Tm- Close to 
but still above Tm, another cluster C{E,V) containing 
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ordered configurations relevant for the crystal state be- 
gins to form, and begins to compete with the disordered 
cluster ViE, V) as far as their entropies are concerned. 
At Tm, there is a transition to the crystal, and the ther- 
modynamics is now controlled by C{E, V). 

As we are interested in supercooled states, we will pre- 
clude the basins that determine C{E, V). This gives rise 
to the restricted PF formalism to study metastability 
[14]. Consequently, we continue to consider 'D{E,V) as 
we reduce the energy. At lower energies or temperatures, 
the cluster 'D{E, V) will be fragmented enough so that 
the percolation ceases to occur mainly because the num- 
ber of components has drastically decreased, and all we 
have are isolated components or clusters of components. 
The disconnected components will then force the system 
to undergo activation-assisted jumps between different 
basins and will force the system to undergo the dynamic 
glass transition. Thus, the picture strongly suggests a 
tantalizing possibility that the above percolation transi- 
tion has direct relevance for the mode-coupling transi- 
tion. We are working on this issue and hope to report 
soon in a separate publication. 



B. Free energy distribution 

The basins that contribute to the thermodynamics are 
enumerated by N'{ip,T,V). All these basins have the 
same equilibrium energy, entropy and the free energy, 
as discussed above. Other basins containing configu- 
rations of energy E{T,V), but not satisfying the other 
conditions do not contribute, even though they may be 
infinitely many, as discussed above. Some of them may 
even have their free energies (p{T, V) <ipb(T, V); they still 
will not determine the thermodynamics because the cor- 
responding free energy F{T, V) due to basins with free 
energy ip(T, V) is higher than the equilibrium free en- 
ergy F{T, V) discussed above in (23). Such a possibility 
{(fij = (f <(ph) has an interesting consequence. The prob- 
ability for probing one of the nonequilibrium basins is 

Vj{T, V) = exp[-5 - - <^b)], (42) 

and is higher than the probability to be in an equilibrium 
basin; here, S is the equilibrium complexity S{iph,T, V). 
There is nothing wrong with it. Let exp[»S(^, T, V)] de- 
note the number of nonequilibrium basins with the basin 
free energy (p. Then, the probability P{T, V) to be in any 
of these basins is obtained by multiplying this number by 
Vj{T,V) above, which yields 

p(r,y) = exp[-/3(i?-F)], 

which surely vanishes for a macroscopic system since F is 
the lowest free energy. Thus, we can neglect the nonequi- 
librium basins in thermodynamic considerations, even 
though they are important for studying transition rates 
from one state to another, which we discuss below. 



C. Transitions among Basins 

We have already argued that all equilibrium basins are 

equally probable; see (34). This is nothing but the Boltz- 
mann hypothesis of microstates [1] ; the only difference is 
that the microstates are now replaced by the basins, and 
the microstate energies E are replaced by the basin free 
energies ip. Since there is no physical interactions between 
basins, each basin can be thought of as representing a 
" microstate" of the system. Thus, the emergence of the 
principle of equiprobability of basins is not surprising. 
What is surprising is that this principle is valid even if 
we consider a restricted PF in which we disallow certain 
configurations (or basins) like those relevant for crystal- 
lization. Such a PF is used to study mctastablc states 
like supercooled liquids [3-5, 14]. Thus, the above princi- 
ple gives rise to a "restricted" ergodic behavior within the 
restricted PF formalism. With this analogy, we can now 
begin to understand the time-evolution of the system, 
since it is well-known [1] that the principle of equiproba- 
bility is equivalent to the principle of detailed balance or 
Liouville theorem: If states (or basins) Bi and B2 have 
the same (free) energy, then the probability 7^12 per unit 
time of a transition from Bi — » B2 and the probability 
V21 per unit time of the reverse transition from B2 — > Bi 
are the same: 

The rate at which the transition Bi — > B2 takes place is 

7^l2 = 7^1 7^12, 

where Pi is the probability that the system explores the 
basin Bi ; we have suppressed the arguments T, V for 
simplicity. It is clear that 7?.i2 = 7?.2i. 

It should be noted that Vij for transitions Bi 
Bj need not all be the same for all choices of i and j. 
Indeed, some of them may be zero. Assume that at a 
given instance, the system is probing a basin Bi and has 
a very small probability of getting into some basin Bq, 
because there is only one pathway through intermediate 
(equilibrium and nonequilibrium) basins B2, B3, B4, .... 
to Bq . While there may be an appreciable probability for 
the transition Bi — ^ B2, Bi B3, etc. because of a 
large number of pathways connecting Bi to B2, Bi to 
B3, etc., the probability of transition Bi to Bq remains 
very small. This also implies that the probability of the 
reverse transition Bq to Bi is also small as there is only 
one path connecting Bq to Bi. It may also happen that 
the probability of transition to other basins from Bq is 
very small. This will mean not only that it is highly un- 
likely for the system to get into the basin Bq, but also 
for it to escape Bq. In no way does this violate the prin- 
ciple of equiprobability. What this situation corresponds 
to is this: once the system begins to probe the basin 
Bq, it stays there for a long time before escaping. The 
amount of time the system probes the basin Bq is still 
the same as the amount of time it spends in any other 
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basins Bi, B2, B^, etc., which is what the principle of 
equiprobability implies [1] . 

The rate Vi at which Vi changes with time is given by 



where we must remember that Vji — Vij . In general, the 
transition probability Vij satisfies 



k^3 



(43) 



The above situation in which Bq has only a few escape 
routes to other basins becomes more probable as the tem- 
perature is reduced. The value of Voj depends, not on 
the free energy heights of the energy barriers but on their 
ratios with T. Thus, at low temperatures, Voj can be very 
small for two disjoint basin components. That does not 
mean that the system will be trapped in Bq, unless there 
are only a few basins (thermodynamically insignificant in 
number) in the sum in (43) so that Vq = 1. li there is 
a large number of basins (thermodynamically significant 
in number) in the sum in (43), T'o < 1 and the system 
will be able to escape to some other basins. It is quite 
clear that the probability of staying in a given basin in- 
creases as T is reduced, mainly because there are only a 
few basins left in the system at that temperature. This 
is the mechanism behind the glass transition. 

To study the trapping of the system to a single basin 
Bq at low temperatures, we can follow the standard idea 
of symm,6try breaking in condensed matter and apply a 
symmetry breaking bias Ao(> 0) to make the free energy 
of Bo lower than the free energies of other equiprobable 
basins: Fq{T, V) = F{T, V)-\q. This will single out the 
basin Bq from all equiprobable basins. At the end of the 
calculation, we can take the limit Aq 0. 



D. Equilibrium identifier Energy £{T,V) 



Let us introduce the mean identifier energy 



£{'P,T,V) = 



Af(ip,T,V) 5 ^ 

' ' j=l,M{ip,T,V) 



£,{V) (44) 



over all basins with the same free energy ip. These need 

not be the equilibrium basins. Furthermore, they not all 
have the same £j{V). Using (42) in (26), we can write 

£{T,V)= J2 £j{V)eM-^-P{'Pj-^b)]. 

j=l,NB{V) 

Let exp[S{ip,T,V)] denote the number of basins with 
the basin free energy ip. Grouping basins according to 
their free energy (p, and introducing F{(p, T,V) = if — 
TS{(p, T, V), we can write 

£{T, V)=Y^ ^(^. T, V) exp[-/?(F - F)]. 



Since F is the lowest free energy corresponding to 'p = 
Pb , we can replace the sum over ip by the most dominant 
term corresponding to ip = (pb, i.e. F = F. This finally 



gives 



£{T,V) = £[pb{T,V),T,V]. 



The sum in (44) for p = ipb is over the equilibrium basins 
counted in J\f{ph,T,V). The equilibrium identifier en- 
ergy must be extensive. The range Tl{T) over which £j 
varies in (44) will depend on the system, and it is con- 
ceivable that TZ{T) /N is not zero in the thermodynamic 
limit, as we discussed above. In this case, the free en- 
ergy fj (T, V) will not be the same for all equilibrium 
basins. However, Goldstein [3] intuitively argues that 
TZ{T)/N should vanish for basin minima as the identi- 
fier, and that the basin minima have a sharp distribution 
about its equilibrium. If this is correct, it would imply 
that all equilibrium basins also have the same average 
minima in additions to the conditions in (S2). This will 
imply that all basins are identical in the thermodynamic 
sense (though not necessarily in the topological sense). 
All we can say is that our analysis is completely oblivi- 
ous of the basin minima or other identifier energies £j. 

Going back to the modified PF in (27), we observe that 
as A — > 0, there is no bias whether lower identifier energies 
are weighted more (A > 0) or higher identifier energies 
(A < 0) are weighted. However, there will in general be 
some topological restrictions. For example, basins with 
Ei^>E{T,V) or E^<E{T,V) are not allowed, but this 
puts no strong restriction on the range TZ(T). For exam- 
ple, all basins with Ei^<E{T,V) regardless of how low 
they are will contribute to E^, (T, V) , as long as the cor- 
responding basin free energy is equal to Ph{T, V). Thus, 
it is not clear whether TZ{T)/N is really zero in the ther- 
modynamic limit. 



E. Complexity as Free Energy? 

The complexity appears in an additive fashion in 
the free energy, which suggests the identification of 
M{p^ T, V) as some kind of a configurational state PF 
ZcsiT,V) [think of Zcs{T,V) as the factor Zbm{T,V) 
in (4)], even though the former truly determines an en- 
tropic quantity S{p,T,V), called complexity; see (7). 
Usually, the microcanonical entropy in statistical me- 
chanics is a natural function of energy, and not temper- 
ature [1], whereas M{p.T, V) is a fimction the tempera- 
ture. Despite this, we now show that it is not proper to 
think of the above identification Zcs{T, V) = M{<p, T, V) 
for the following reason. The identification will yield 
the free energy Fcs{T,V) = -T\nZcsiT,V) = -TS, 
from which we can compute the corresponding entropy 
[-{OFcsiT, V)/dT)v], and energy [E = F + TS] 



S + T (dS/dT^^, 



Scs{T,V) 
Ecs{T,V)=T^ (dSjdT^ 
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We now use this in (23), and obtain 

F{T,V) = E^{T,V) + Ecs{T,V) 

- T[Sb{T, V) + S{T, V) + T (^dSjdr'^ 

This form of the free energy will give for the energy the 
first two terms, which is not correct; see (39). Similarly, 
the quantity in the square rackets also is not the correct 
entropy of the system; see (7). Hence, complexity should 
not be considered as a free energy. Its proper identifica- 
tion is with the entropy. This justifies using the complex- 
ity as the entropy associated with configurational states. 
This entropy is neither the entropy Sqs {T, V) nor the en- 
tropy Sbi{T, V) due to any basin identifier. The entropy 
associated with the configurational states advocated by 
Goldstein truly represents the entropy associated with 
the basins of a given free energy (/?, and not a given BM 
energy. Moreover, the configurational states of Goldstein 
should be identified as the basins of a given free energy 
ip. 

F. Extension to NTP-Ensemble 

For the TVTP-ensemble, we need to integrate 
Z{T,V)exp{-pPV)dV, see the PF in (1), over the vol- 
ume, so that the role of V in the above discussion is 
played by the equilibrium volume V(T, P) and should 
be replaced by it. All the above formulas remain valid 
if wc replace V by either P or y(T, P) in all thermody- 
namic quantities. (It is for this reason that wc have ex- 
plicitly shown V in all the formulas above.) For example, 
the basin energy is represented as a function Ej (T, P) or 
Ej [r, V{T, P)]. The same is true of the basin free energy 
(fi, which now represents the Gibbs free energy. The com- 
plexity is to be treated as a function of ip, T, and P or 
F(T,P). 

The PE landscape is now defined in a (3A^ + 2)- 
dimensional landscape space Ajv in which V is added 
as an additional axis. The PE surface is now a (3A'' -|- 1)- 
dimensional surface E. Since the configuration space 
T[V) requires a particular volume V for its definition, we 
need to project S for a particular value of V. Thus, we 
need to deal with the same sort of projection n(V) in the 
configuration space as before. The slice dT, also needs to 
be projected for a particular value of V. Thus, we deter- 



mines the number of configurations W{E, V) for a given 
E, and V. In equilibrium, corresponding to E, and V exist 
T, and P such that E{T, P) = E, and V{T, P) = V. All 
the conclusions above remain valid. 



G. Summary 

In summary, we have investigated the landscape de- 
scription at all temperatures, even though the picture 
was originally introduced by Goldstein as valid and use- 
ful only at low tempertures. We have succeded in val- 
idating the picture at any temperature. Thus, the re- 
sults we present here are not restricted to any particular 
temperature. We have shown that the basins represent 
the configurational states that were introduced by Gold- 
stein. These basins and their free energy (p play the role 
of microstates and their energy E in the microcanonical 
ensemble, and are governed by a similar equiprobabil- 
ity principle, which extends the Boltzmann equiprobabil- 
ity principle for microstates. The corresponding entropy 
due to the number of basins of a given free energyy is 
called the complexity [9], and correctly represents an en- 
tropy and not a free energy. We have shown that the 
entropy is a sum of complexity and intrabasin entropy. 
We also conclude that all equilibrium basins at a given 
temperature not only have the same equilibrium free en- 
ergy i^b(r, y), but also have the same energy E{T,V), 
and the same intrabasin entropy S\^{T,V). Because of 
the relevance of equilibrium basin identifiers in different 
temperature ranges,we have provided a procedure to ex- 
press thermodynamic quantities as a function of equilib- 
rium basin identifier energy £{T, V) since the equilibrium 
free energy does not explicitly depend on £{T,V). Our 
analysis provides justifications for two of the conjectures 
made by Goldstein but with slight modifications. The 
equiprobability theorem allows to draw some conjectural 
conclusions about the dynamics and the role of landscape 
projection components in the phase space. We discuss 
the percolation of basin projection components and spec- 
ulate about the possible connection of the disappearance 
of percolation with the dynamic glass transition. This 
topic is being investigated at present and we hope to re- 
port on it in near future. 
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